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I. INTRODUCTION 



Much progress has been made in the study of low dimensional systems of strongly cor- 
related electrons. A notable example is the one dimensional Hubbard model, which was 
solved by Lieb and Wu [I| using the Bethe-ansatz. This model exhibits a hidden SO (4) 
symmetry m which is the basis of its properties. The suggestion of Anderson that the 
two dimensional Hubbard model contains the essential physics of high temperature super- 
conductivity and that its normal state may share the Luttinger-liquid-like feature of one 
dimensional interacting electron systems ||,|5| has induced renewed activity in the study of 
such strongly correlated one dimensional models. 

Recently, an extended Hubbard model with nearest neighbor interactions was introduced 
IK?], and is solvable by Bethe-ansatz in one dimension ||. The spin SU(2) along with the 
^-pairing SU (2) of this model are combined and extended to form a SU (2\2) supersymmetry. 



We have introduced the Calogero-Sutherland @,[T(J version of the extended Hubbard model, 
defined both on a uniform lattice and on a non-uniform lattice in one dimension. The long 
range versions of the extended Hubbard model are exactly solvable, and the ground states, 
excitations and integrabilities have been studied in considerable detail |TT|JT2] . 



In this paper, we generalize the previous results by extending the supersymmetry to 
SU(m\n). Thus we introduce the SU(m\n) strongly correlated electron system of 1/r 2 long 
range interaction. This general model encompasses the SU(0\2) Haldane-Shastry spin chain 
13| , [T4l| and the S77(l|2) supersymmetric t-J model fl5| as well as the SU{2\2) extended 



Hubbard model. 

In the next section we present the general features of the SU(m\n) supersymmetric 
model. Following this, we study the system, first defined on a non-uniform lattice given 
by the roots of the Hermite polynomials, and then defined on a uniform lattice. In both 
cases, the SU(m\n) models are completely solvable. We explicitly construct the ground 
state wavefunctions and provide the ground state energies. We also examine the excitation 
spectra of these systems, making use of their integrability properties. 
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II. THE SU(M\N) SUPERSYMMETRIC ELECTRONIC MODEL 



Consider a system defined on a lattice with a total of L sites. At each site there is exactly 
one particle, which may either be a fermion with n flavors (internal degrees of freedom) or 
a boson with m flavors. Since there are a total of m + n states at each site, n fermionic and 
m bosonic, this is the basis of the SU(m\n) supersymmetry. 

In order to describe these degrees of freedom, we use a representation in terms of fermion 
and boson operators, / and b, satisfying the usual (anti-)commutation relations 

{fia, fja'} = dijtiaa' {fi<r, fja'} = 

ha, b] a ,] = 5ij5 aa > [bia, b ja >\ = . (1) 

The indices i and j label the site, a — 1, . . . , n labels the fermion flavor and a = 1, . . . , m 
labels the boson flavor. These operators may be represented in a supersymmetric manner 
by c iv = {fia}, {ha} where v — 1, . . . , m + n labels the (either fermionic or bosonic) species. 
The c's then satisfy the ( ant i-) commutation relation 

[c iu ,c] u ,]± = 5ij5 vu/ , (2) 

where [. . .]± indicates either a commutator or an anti-commutator as appropriate. In general, 
this representation allows for multiple occupancy at each site. In order to project onto single 
occupancy we use the Gutzwiller operator, Pq, defined by 

P G = U^E etc 

1 = 1 

L 

= II 5 iy fi /. +y &t b. ■ ( 3 ) 

1 = 1 

The Hamiltonian for the SU(m\n) supersymmetric model takes the form 

H = ~\pg E Jid^aPa, (4) 

where the coupling parameter Jij = l/d(i,j) 2 takes an inverse square form. The distance 
function, d(i,j), depends on the lattice and is what distinguishes between the models on 
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uniform and non-uniform lattices. In this section, we concentrate on the features of H 
that are independent of the specifics of d(i,j). The key features of this model, especially 
its integrability, are derived from the graded permutation operator, IT,,,-, which exchanges 
particle species at sites i and j. In the above representation, the permutation operator takes 
the form 

m+n m+n 

v=\ v'=\ 

n n mm 

= ~ £ £ fLfia' fja'fja + £ £ ^ia^ia'^ja'^ja 
cr=l <t'=1 a=l a'=l 

n m 

+EEW<4+M • ( 5 ) 



Since the Hamiltonian only permutes particles, the number of bosons and fermions of 
each flavor are all conserved quantities. It is also easy to verify that H commutes with both 
the SU(n) fermion and SU(m) boson flavor symmetries generated by 

S a = "t S?, S? = £ flT a aa , {n)f i(T , , (6) 

i=l <t,<t' 

and 

V a = J2 Vt, V? = E *>lT a aa > {m)b ial , (7) 

i=l a,a' 

respectively (here T(n) and T(m) are hermitian SU(n) and SU(m) generators in the fun- 
damental representations). In fact, ignoring the grading, since Ely acts symmetrically on 
all m + n species [as is evident from the first line of (Bp], we would expect a SU(m + n) 
symmetry of this model. Taking statistics into account, this is modified into a SU(m\n) 
supersymmetry of which the bosonic subgroup is given by SU(m) x SU(n) x U(l) (see e.g. 
Ref. [0). The U(l) symmetry is related to the boson and fermion number currents, and is 



given by 

L 



Lj -i n -r m 

j = Y,^ 4 = -£/i/* + -£ 6 U a - (8) 



i=l cr=l a=l 



Although it appears there is another U(l) symmetry given by the remaining orthogonal 
combination of boson and fermion currents, this is actually not the case since we restrict 



ourselves to single occupancy, and thus total particle number is fixed. This means that the 
second U(l) symmetry is trivial (except in the special case when m = n). 

Complementing the above bosonic generators of SU(m\n) are the fermionic generators 

L 

Qua ^ t Qiaay Qiaa fia^ia 1 

i=l 
L 

Qaa Qiaui Qiaa ^ictficr j (9) 

i=l 

transforming as the (m,n) and (m,n) of SU(m) x SU(n) respectively. The complete 
SU(m\n) superalgebra is given by these 2mn fermionic generators as well as the m 2 + n 2 — 1 
bosonic generators of Eqns. (0), ([?P and 

In general, since the Hamiltonian H is SU(m\n) invariant, all states in the system fall 
into complete representations of the SU(m\n) supersymmetry algebra. These supermultiplets 
have a very rich structure, and span states with different numbers of fermions and bosons. 
Representations of SU(m\n) may be divided into two categories, "typical" and "atypical" 
|T6| . The typical representations group together states in sectors with boson and fermion 
numbers given by (Q, M) = (6, /), {b + 1, / — 1), {b + 2, / — 2), . . . , {b + mn, f — mn), while 
the atypical ones terminate early and hence span fewer than the mn + 1 sectors of typical 
ones. In both full SU(m\n) representation may be decomposed in terms of a 

direct sum (over sectors of decreasing fermion number) of representations of the bosonic 
subgroup SU(m) x SU(n) x U(l) C SU(m\n). For an irreducible representation, we refer 
to the state with greatest fermion number (i.e. the (b, f) sector) as the highest member of 
a supermultiplet and the state with the fewest fermions (the (b + mn, f — mn) sector for 
typical representations) as the lowest member. 

While the important properties of the system may be understood through this super- 
multiplet structure, we generally wish to have greater control over the individual species 
populating the lattice. In particular, we note that since SU(m\n) has a maximal Cartan 
subalgebra given by [/(l) m+n_1 , there are a set of m + n — 1 conserved currents which are 
equivalent to the fermion and boson number currents of the m + n species (along with the 
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single occupancy constraint). As a result, the number of bosons of each flavor and the num- 
ber of the fermions of each flavor are all conserved quantities, and hence we may work in a 
subspace where all these quantities are fixed. This is easily accomplished by adding to H a 
set of chemical potentials for the different species: 

m+n L 

H = H+J2» {l/) T, c l c i^ ( 10 ) 

V=l 1=1 

(one of the /z's is redundant). Each subspace of fixed occupancies may be labeled by the 
m + n numbers {Q a , M a } given by 

L 

M a = J2fiJi*> <r = l,2,...,n 
i=i 

L 

Qa = J2 b Lb ia , a = l,2,---,m . (11) 

i=l 

Since each site is singly occupied, these quantities satisfy the relation that J2a M a + J2 a Qa = 
L. Working in such a subspace, we have broken the SU(m\n) symmetry down to its maximal 
Cartan subalgebra. Alternatively, we often consider the case when SU{m\n) is not broken 
completely down to C/(l) m+n_1 , but instead to its bosonic subgroup SU(m) x SU(n) x U(l). 
This is equivalent to working in a subspace of fixed (Q,M) where Q = J2 a Qa is the total 
number of bosons and M = J2a M c is the total number of fermions, and corresponds to 
using the U(l) current J to break supersymmetry. 

We now define the wavefunctions in this Hilbert subspace of fixed ({Q a }, {M a }) by 
writing the state vectors of the system in the following form: 

|0)= £ Hxa;y a )UfUIlKua^) ■ ( 12 ) 

{xa},{ya} 1=1 k=l 

The amplitude (f)(xa, ya) is antisymmetric in the fermion positions and spins, {xcr} = 
(xi<7i, • • • , Xm&m), while symmetric in the boson positions and spins, {yea} = 
(yicxi, ■ • ■ , Vqoiq). Because of the supersymmetry, it is convenient to combine the boson 
and fermion positions by defining (gi, q 2 , ■ ■ ■ , q£) = (x±, • • • , Xm\vi, • • • , Uq) and (f)(xa, ya) = 
<f)({q}, {&}, {a})- Due to single occupancy, the set {q} spans the lattice. 
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With the state vector written in the above manner, the fermionic and bosonic nature of 
the particles is encoded in the symmetry properties of the wavefunction. Thus the graded 
permutation operator, (|5]), takes a particularly simple form independent of the particle 
statistics. The resulting Hamiltonian, acting on wavefunctions 4>({q}, {cr}, {a}), is 

= "2 1 <fe L %^F' (13) 
where My interchanges particles i and j, Mjj0({g}, {a}, {a}) = <fi({q'} , {a} , {a}) , with 
{<?'} = (Qi, Q2, 4 4 ■ , Ql) = (<?i, • • • ,Qjr ■ ■ ■ ■ ,Ql)- The resulting eigenenergy equation then 
takes the form 

-J E d(q i ,q J )- 2 M lJ <j ) ({q},{a},{a}) = E<j ) ({q},{cr},{a}) . (14) 

We immediately see that the lowest energy state in the full Hilbert space corresponds 
to constant 0. Since this is only compatible with bosonic symmetry, we conclude that the 
SU{m\n) ground state is in the M = sector and is described by the wavefunction 

<t> G {ya) = 1 . (15) 

Because this is a good wavefunction independent of boson flavors a, the ground state is in 
general degenerate, and is given by the L-fold symmetric combination of the fundamental 
representation of SU(m). In terms of the Young tableau, this is the ([L 1 ], [•]) of SU(m) x 
SU(n). Since the wavefunction is constant, the ground state energy is given by 

E G = -\ E d(i,j)-\ (16) 

and may easily be summed. All excitations, whether fermionic or bosonic, are built on top 
of this. 

To be more precise about the nature of the ground state, the purely bosonic (Q, M) = 
(L, 0) state is just the lowest member of a SU{m\n) supermultiplet, and therefore its 
degeneracy in the full Hilbert space is larger than indicated by the M = sector 
alone. We find that the ground state representation is atypical, and spans the sectors 



(L — n, n), (L — n + 1, n — 1), . . . ,(L,0) with corresponding representations ([(L — n) 1 ], [l n ]), 
{[{L— n+1) 1 ], [l n_1 ]), . . . , ([L l ], [■]). Alternatively, this result follows by noting that fermionic 
wavefunctions which are symmetric in the coordinates {x} may be constructed by taking 
antisymmetric combinations of their flavors a. 

The highest energy state of the SU(m\n) model is given by a completely antisymmetric 
wavefunction, corresponding to the fully fermionic sector, M = L. In the absence of periodic 
boundary conditions, the wavefunction is 

(j) M (x(j) = Y[(qi - qj) , (17) 

i<j 

with energy Em = —Eg- Since <pM is independent of the fermion flavor a, the highest 
energy state corresponds to the ([•], [L 1 ]) representation (this time the highest member of a 
supermultiplet). As a result, all states in the spectrum lie in the energy range [E G , —Eg]- 
Furthermore, any eigenstate of H, when multiplied by 4>m, remains a good eigenstate, but 
with opposite eigenenergy. Since this also interchanges bosonic and fermionic boundary con- 
ditions, it corresponds to the interchange of the SU(m\n) and SU(n\m) theories. Therefore 
this indicates that H for the SU(m\n) model corresponds to —H for the SU(n\m) model 
T7J. We will subsequently make use of this symmetry to obtain the upper bound of the 



energy levels of the SU(m\n) permutation model defined on the Hermite lattice. 

In general, we wish to work, not in the full Hilbert space, but rather in a given subspace 
of fixed (Q,M). In this case, wavefunctions have mixed symmetry properties. We will 
consider the specifics of these wavefunctions in the next two sections. Here we mention 
that, since the ground state in this sector contains a symmetric combination of Q bosons, 



it is again degenerate [0, and corresponds to the Q-fold symmetric combination of the 
m of SU(m). As a general rule, in order to lower the energy, the wavefunction should be 
as symmetric as possible. This may be accomplished by distributing the M fermions as 
equally as possible into the n flavors, resulting in the ground state representation ([Q 1 ], [V]) 
where r = Mmodn. This will be made clear subsequently for both the non-uniform and 
the uniform lattice. 
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III. SU(M\N) ELECTRONIC MODEL ON A NON-UNIFORM LATTICE 



As we have seen above, many properties of the SU(m\n) model are quite general, and 
result from the permutation nature of the system. In order to proceed further with a detailed 
study of the excitation spectrum, we now fix the lattice, as specified by the distance function 
d(i,j). In this section we discuss the SU(m\n) model defined on a non-uniform lattice in 
one dimension where the sites of the chain are given by the roots of the Hermite polynomial 
Hl(x). It is well known that this Hermite polynomial has L roots, r 1; r 2 , . . . , r^, which are 
all real and distinct. Thus the distance function for this non-uniform lattice is defined by 

d(i,j) = \n -Tj\ . (18) 

This one dimensional chain is well defined and yields an exactly integrable system. The 
SU(m\n) model on a non-uniform lattice encompasses the SU(0\n) spin chain with inverse 
square exchange |L{| and the SU(l\n) supersymmetric t-J model as well as the 



SU{2\2) extended Hubbard model JTT] 



Since the Hamiltonian, (|13|), is written solely in terms of permutations on the particle 
positions and is independent of statistics, many results for the SU(0\2) spin chain on a non- 
uniform lattice |19| are applicable to the general SU(m\n) system as well. In particular, the 
Hamiltonian commutes with a set of conserved quantities {I a}, namely 

[H,I A ] = 0, [I A ,I B ]=0, ,4,5 = 0,1, 2,... ,oo , (19) 

where 

i=l 

L 

4 = * (9k- qjY l M kj + iq k , a k = (4) f . (20) 

These relations yield the integrability conditions of the SU(m\n) model. The quantities I a as 
written act only on the coordinates of wavefunctions (fi({q}, {&}, {a}) and thus commute with 
the SU(m\n) algebra. One can use the permutation symmetry property of the wavefunction 
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defined by Eqn. (|Ll|) in order to write explicit representations for I a in terms of the b and 
/ operators. 

As shown in the previous section, the ground state of this model in the complete Hilbert 
space has its lowest component in the (Q, M) = (L, 0) subspace, and has an energy 

^-Ui^'-i 1 ' 1 - 1 '' (21) 

While this is the ground state of the system in the complete Hilbert space, we now wish to 
investigate the model in a given subspace of fixed occupation numbers ({Q a }, {M a }). Based 
on our previous experience with the supersymmetric t-J model and the extended Hubbard 
model, we anticipate the ground state wavefunction in this sector to take the following form: 

(f) (xa,ya) = J] ~ xj) 5 ^ e^ 11 ^"^ . (22) 

l<i<j<M 

This is essentially the minimal possible wavefunction that satisfies the appropriate anti- 



symmetries under fermion exchange. Using the techniques of fl20,21], we may prove that 



this Jastrow wavefunction is an eigenstate of the Hamiltonian ( |T3"D with eigenenergy 

1 1 n 

E = —L{L - 1) + 2 E M °( M ° ~ 1) > (23) 

<T=1 

independent of bosonic species, as anticipated in the previous section. Minimizing this 
energy subject to the constraint Yj<t M a = M gives 

E(Q, M) = --ML - 1) + ^[M(M - n) + r{n - r)\ , (24) 
4 In 

where r = M mod n. The ground state wavefunction has the fermions distributed as evenly 
as possible among the n distinct flavors. For M > n, this is the highest component of the 
supermultiplet. 

We now turn to the question of excitations above the ground state. As in the super- 
symmetric t-J model and the supersymmetric extended Hubbard model, one can show that 
there are several ways to create excitations from this ground state. The first way is to 
excite the M fermions, the second way is to excite the Q bosons. In fact, a k and a k are 
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raising and lowering operators of the Hamiltonian H [pO| , pT| . To construct excited state 
wavef unctions, we first use the convention implicit in (^) that the fermions are ordered 
before the bosons. Thus the positions of the fermions are 52, • " ' > Qm) and the positions 
of the bosons are (qM+i, • • • , Ql)- One way to excite those / fermions is to construct the 
following wavef unctions: 

_ M 

\K t , K 2 ,...,K M )= Y, Ii( a l°i) KPi \<h) , (25) 

Pi,P 2 ,--,P M i=l 

where the summation P is over all possible permutations, Ki are nonzero integers, and |0o) is 
the state vector corresponding to the amplitude of Eqn. (|22). These states, if not vanishing, 



will be eigenstates of the Hamiltonian, with eigenenergies given by 

M 

E{K 1 ,---,K M )=E + Y K i ■ ( 26 ) 

i=i 

To excite the bosons, one way is to construct the state vectors: 

Q 

\K U K ar ~,K Q ) = Y, Il( a I + M^)^ I0o) , (27) 
p i=i 

with eigenenergies given by 

Q 

E = E + Y K i 1 (28) 

i=l 

if the state vectors constructed this way do not vanish. There are also some other ways to 
excite the fermions and bosons, such as 

\K h K 2 ,---,K Q )=Y1 flKU)^ I0o) • (29) 
p i=i 

Now, let us return to the question of why <p is anticipated to be the ground state of the 
system. Heuristically, we note that, by construction, <f) is as symmetric as possible while 
still compatible with fermion statistics. Thus it is as close as possible to the maximally 
symmetric state 4>g- More rigorously, we can prove explicitly that 

Q Q 

(Y a i+M(Xi)\(i>a) = 0, (Y ««+m) |0o) = . (30) 
i=i i=i 
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Furthermore, one can show that 

M M 

(E a ^)l<M = 0, (E a Ol0o) = O, (31) 

i=l i=l 

for the fermions' degrees of freedom. These relations yield the impossibility of constructing 
non- vanishing states of energy E — 1 by using the lowering operators a^. We therefore can 
regard these identities as partial confirmation that 0o is the ground state. A more complete 
proof requires further work. 

In the subspace of fixed ({Q a }, the full energy spectrum of the SU(m\n) model 

on this non-uniform lattice is expected to consist of equally-spaced energy levels: 

E s = E + S , (32) 

where S = 0, 1, 2, • • • , S max . There is an upper bound on the values of S due to the finite 
size of the Hilbert space. Using the boson-fermion interchange symmetry, H — > —H for 
SU(m\n) — > SU(n\m), and Eqn. (f?3|), we find that E Smax = \L(L - 1) - \ T, a Qa{Qa - 
1). Presently we are unable to characterize the pattern of the energy spectrum by some 
systematic rule, nor are we able to explain the degeneracies of the energy levels by the 
underlying symmetries of the system (presumably Yangian symmetry). 



IV. SU(M\N) ELECTRONIC MODEL ON A UNIFORM LATTICE 

We now turn to the SU(m\n) supersymmetric electronic models on a uniform lattice in 
one dimension. This uniform lattice is characterized by the trigonometric interaction 

d(i,j) = £smf?^n ■ (33) 

The general SU(m\n) model reduces to the Haldane-Shastry spin chain p~3| , |T4fl for the case of 
SU(0\2), the supersymmetric t-J model introduced by Kuramoto and Yokoyama p!5| , |22|j23[] 
for SU (1\2), and the extended Hubbard model studied by us previously |12|] for SU(2\2). 
For the Hamiltonian ( |T3"D on a uniform lattice, we may use the results of Fowler and 



Minahan for the Haldane-Shastry spin chain |24| , except that one has to use the permutation 
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properties of the amplitude as defined by Eqn. fll2|) . The result is a set of conserved 
quantities I A , 

[H,I A ] = 0, [I A ,I B ]=0, A,B = 0,l,2,...,oo, (34) 

where 

1=1 

^ = £ / y Z % 7 %. ^ = ex P(^r^) > ( 35 ) 
fc(^-)=i ~ z j'J ^ 

which yield the integrability of the SU(m\n) electronic model. 

In the full Hilbert space, the ground state (as the lowest component of a supermultiplet) 



is obtained when M = 0, with trivial ground state wavefunction (15[). The corresponding 
ground state energy on the uniform lattice is 

E G /^/Lf = -\ J2 ■ 2 , r ^7n = -\l(L 2 ~ 1) , (36) 

with degeneracy of the lowest component given by the ([L 1 ], [•]) representation of SU(m) x 
S?7(n). 

Once again, we turn to the examination of the ground state in a given subsector spec- 
ified by ({Q a }, {M a }). As before, when fermions are present, we expect the ground state 
wavefunction to have the minimum antisymmetry possible consistent with fermi statistics. 
This time, however, due to the uniform lattice, we also require the wavefunctions to have the 
appropriate periodicity under q { — > q^ + L. The results for the SU(l\n) case |25j motivate 
us to write the following Jastrow wavefunctions for the SU(m\n) case: 

M 

o (xa, ya) = J] xf'* J[ {X t - X s )*°*'i C W^) , (37) 

i=l l<i<j<M 

where = exp( 2 ™ fc ), and M = Y,a M a is the total number of fermions. These wave- 
functions are eigenstates of the Hamiltonian provided the quantum numbers Ji, J2, • • • , J n 
(uniform lattice momenta for each fermion flavor) obey the constraints — M a < J a < 0. 
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The eigenenergies corresponding to the above wavefunctions were derived in |[25|| , and 
take on a more concise form when expressed as a function of the shifted momenta, K a = 
J a + h(M a — 1). Without loss of generality, we assume the occupation numbers, {M a }, are 
ordered according to Mi > M2 > ■ ■ • > M n . In this case, we find 

11™ 1 
E = --L(L 2 - 1) + - £(3L + (a - 2)M a )(M 2 -1)--J] M a {M 2 - 1) 

<J 1 CT<cr' 

n 

+2(L - M) £ ^ + 2 X; - i^) 2 , (38) 

(7=1 (7<<T' 

provided the (either integer or half integer) K^s satisfy 

\K.\ < - 1) 
\K a -K a ,\< l -(M a - M a ,) for a 1 > a . (39) 

The lattice momentum carried by this state (over the bosonic background) is simply Pq = 

Examination of Eqn. ( |38|) indicates that the lowest energy state in the above set of 
Jastrow wavefunctions is reached when the K^s are all as close to as possible. This 
is accomplished by taking K a = for odd M a or 1/2 for even M a . Note that whenever 
some of the M CT 's are even, we could equally well have taken K a = —1/2, leading to a 
two-fold degeneracy of the ground state arising from the reflection symmetry of the lattice 
(in addition to the degeneracy arising from the bosonic species). As in the case of the non- 
uniform lattice, the ground state in a given (Q, M) sector has the M fermions distributed 
as evenly as possible. The ground state energy is then given by 

E(Q, M)/(n/L) 2 = -\h[L 2 - 1) + -L(3L - M)(M 2 - n 2 ) 

on 

+ ^-r(n - r)[3(L + M) + 2(n - 2r)\ 

1 ( n — r, if Int (M/n) is even; 
+-(L-M)\ (40) 

[ r, otherwise. 

As before, r is given by r = M mod n. Only the first line is important in the thermodynamic 
limit. This ground state (and its reflected pair for P ^ or L/2) transforms as the {[Q l ], [l r ]) 
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of SU(m) x SU(n) (once again the highest component of a supermultiplet for M > n). We 
note that when all M a are odd, this particular wavefunction can also be obtained by taking 
the strong interaction limit of the corresponding continuum quantum system | [26|| . 

The full energy spectrum for this long range permutation model may be derived using the 
asymptotic Bethe-ansatz (ABA) [|10|| . The motivation for this ABA lies in the fact that the 
scattering is essentially two-body in nature, even in the presence of long range interactions, 
as indicated by the Jastrow form of the wavefunction. The ABA has been used successfully 



in the S77(l|2) supersymmetric t-J model p7 as well as its SU(l\n) generalization |[22|| , and 



gives exact results, even in the non-asymptotic regime, as proven to be true for the SU(0\2) 
Haldane-Shastry spin chain |28) and for the SU(1\2) t-J model [^3| — a fact that is related 



to the integrability of the system. 

For the SU(m\n) generalization of the ABA, we need to treat multiple flavors of both 
bosons and fermions on the lattice. In the case of nearest-neighbor interactions, Lai first 



discussed such lattice permutation models with mixtures of bosons and fermions in 1974 [29 



This was later generalized by various authors [|T^,^,^,g. For the multi-species case, the 
Bethe-ansatz takes on a nested form, where one species is removed at each step. For the long 
range SU (m\n) model, this is accomplished by starting with one of the m + n (either bosonic 
or fermionic) species as the background, and then peeling off the remaining m + n — 1 species 
one at a time, resulting in m + n — 1 sets of Bethe-ansatz equations. There are a total of 
(m + n)l ways to perform this nesting. However this number may be reduced by the obvious 
nl fermion permutations and ml boson permutations, yielding (m + n)l / (mini) independent 
ways of writing the nested Bethe-ansatz. A particular choice of nesting may be denoted 
by a string of m + n 5's or F's, indicating the order in which the bosonic and fermionic 
species are removed, working from right to left (the rightmost character denotes the choice 
of background). Of course all such choices of the nesting should yield equivalent results. 
The ABA obtained previously for the SU(l\n) supersymmetric t-J model |22| corresponds 



to the BF n choice of nesting. 

For a general v (= m+n) component system, we encode the nesting of the ABA according 
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to the statistics Sj of species % where s« = for bosons and 1 for fermions. Since we take the 
first species as the background, the B m F n nesting corresponds to {s} = (1, . . . , 1, 0, . . . , 0). 
We denote the occupancies of the species by Mj (here i runs from 1 to v) and define the 
quantities 



3 i 



(41) 



j>i 



where Nq = L and N u = 0. For the nested ABA, we introduce v — 1 sets of pseudomomenta, 
{p^}, {p^}, • • • , {p^ 1 '} where each set {p^} consists of the N a quantities 



P 



(a) 



i = 1,2 



, ii, . . . , i i a , 



(42) 



all of which are within the range [— it, n]. The energy and the total lattice momentum of the 
system depend only on the p^'s, and are given by 



E 



P 



\Sl+l 



7T 



1 , 1 



B *(i-fi) + §E«*! 



i=l 



(i) 



7T 



8=1 



mod 2ir . 



(43) 



The pseudomomenta, {p^} for a = 1, . . . , v — 1, are obtained by solving the v — 1 set of 
equations of the nested ABA, expressed concisely in the following form: 

iVa + l 



N a -! N a 

E %S a) - pt 1] ) = 27r/| a) + 5 Sa , Sa+1 E %S a) - 
fe=i j=i 



i=i 



(44) 



where 9{x) = 7rsgn(a;) is a step function. For a = 1 the left hand side of (0) is replaced by 
Pi^L, and for a = v — 1 the last term on the right hand side is dropped. A given state in the 
ABA spectrum is thus specified by the set of quantum numbers {1^}. At a given level a of 
the nesting, there are N a non-overlapping /■"'' s (all integers or half-integers as appropriate) 
which are required to lie in the range 



l4 a) l < ^(^Vi-^a+^i-i) 



(45) 
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for s a = s a+1 or 

l4 a) |<^(A^i-iV a+1 -2), (46) 

for s a 7^ s a+ i. These conditions lead to the restrictions on the occupation numbers, either 
M a > M a+ i or M a > N a+1 + 1 respectively. An exception happens when N a = 0, in which 
case the nested ABA terminates early, giving no further restriction on M a . 

We remark that in general the ABA still provides exact eigenenergies when conditions 
([4"5|) and fllB]) are relaxed by one unit on either end. However we have found the above 
restrictions to be necessary in order to ensure that the ABA generates only highest weight 
states of the SU(m\n) superalgebra. It is important to realize that the choice of nesting 
plays a crucial role in defining what is meant by the highest weight representation. In 
particular, for SU(m\n), the Cartan subalgebra lies in the bosonic subgroup SU(m) x 
SU(n) x (7(1) and has rank m + n — 1. Thus highest weight representations of SU(m\n) are 
annihilated by a set of m + n — 1 simple positive roots. It is this set of simple roots which 
is determined by the nesting of the ABA, with identical statistics in the nesting (s a = s a+ i) 
corresponding to bosonic roots and s a ^ s a+ i corresponding to fermionic roots. For this 
reason the most straightforward nesting is given by B m F n (or F n B m ), in which case there 
are m + n — 2 bosonic roots corresponding to the simple positive roots of SU(m) x SU(n) and 
a single fermionic root taking a state from the (Q, M) sector to (Q — 1, M + 1). This gives 
the conventional definition of the highest weight state of SU(m\n) residing in the highest 
component of the supermultiplet. 

In order to be more specific, we now pick the natural B m F n nesting and the original 
notation of M a and Q a for fermionic and bosonic occupation numbers respectively. The 
nested ABA may be written in terms of the n — 1 sets of "fermionic" pseudomomenta 

n 

pt ] : t = l,2,...,N a where N a = Q + £ M a , (47) 

o=a+l 

and m sets of "bosonic" pseudomomenta 

m 

qf: fx = l,2,...,N' b where = £ Qa • (48) 

a=b 
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The ABA may now be written explicitly as 

P?L = + £ ^ ~pP) ~ E^ -Pf) , 

E } - P?>) = 2vr/f + E } - pf ) " E ) - Pf ) 



Y,^t l) -p { r 2) )=^it l) + e o(p i r 1) -ptr 1) )-Y,o(pt 1) -$ ) ) 
j fe'(^fe) ^ 

E %? } - Pt 1] ) = 27T jw + E } - ^ 2) ) , 



1/i 

fc 



E %1 2) - } ) = 2-4 2) + E %1 2) - i 2) ) - E %1 2) - ?r 



E%7 m_1) -^ 2) ) = 2 ^ 1) + E %f- 1) -^~VE%f~ 1) -^ ) : 

7'(^7) 1 



E %S m) - ^ m ' 1} ) = 27r 4 m) + E %S m) - ^ ra> ) - (49) 

7 r/^r?) 

with quantum numbers {/( a )} and {J^}. It is anticipated that the ABA span the full 
set of energy levels when both the highest weight property of the ABA states and the 
supermultiplet structure connecting states in different (Q, M) sectors are taken into account. 
We have verified that this is the case numerically for the SU (2\2) model on small lattices. The 
numerical results also indicate that, while the ABA accounts for the full energy spectrum, it 
does not give all the proper degeneracies of the states. The missing states presumably arise 
from localized bound states which are not described by the asymptotic Bethe-ansatz. Thus, 
in order to properly count the degeneracies and derive the full SU(m\n) supermultiplet 
structure of the spectrum, one must generalize the "squeezed string" picture of the SU(n) 



Haldane-Shastry spin chain ||32 ] . 



V. CONCLUSION 



In this work, we have introduced the SU(m\n) electronic model of long range interaction, 
which is integrable on both uniform and non-uniform one dimensional lattices. The ground 
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state, excitation spectrum, and integrability properties of this model have been considered 
in detail. Since the SU(m\n) model includes many previously well-known models as special 
cases, we have shown how a uniform treatment of all such models may be accommodated. 

For the integrable SU(m\n) model on a non-uniform lattice, perhaps the most interesting 
open problem is to explain the high degeneracy of the equally-spaced energy levels using 
the underlying symmetries of the system. Presently it is still unclear how to construct the 
thermodynamics of the long range SU(m\n) permutation model on the Hermite lattice. 

For the uniform lattice, since the energies of the ABA states are given by the essentially 
non- interacting formula, (|43|) , we note the interesting result that for a given lattice of length 
L, the set of allowed energy levels are completely determined, and is independent of the 
number of fermionic and bosonic species, n and m [jlT]] . The differences between the SU(m\n) 
models on a uniform lattice hence lie only in the degeneracies and supermultiplet structures 
of the models, with larger n and m leading to higher degeneracies. 

Since the full SU(m\n) symmetry is broken in most physical problems, we generally wish 
to work in a given subspace of fixed (Q, M), corresponding to the breaking of supersymmetry. 
It is clear that, in this case, a proper understanding of both the supermultiplet structure 
(including degeneracies) and SU(m\n) representation theory is required. We note that the 
latter is quite intricate as atypical representations are clearly involved. 

Working with fixed (Q, M), one only sees the bosonic SU{m) x SU(n) subgroup of the full 
SU(m\n) symmetry group. For both the non-uniform and uniform lattice, we have explicitly 
found the ground state wavefunctions and their corresponding eigenenergies. Because of the 
symmetry properties of the wavefunction, the ground state wavefunction is independent of 
the bosons. Hence the ground state in any Q > sector is degenerate whenever m > 1, and 
transforms in the ([Q 1 ], [l r ]) representation where r = Mmodn. 

The ground state structure and ABA results for the excited spectrum of states on the 
uniform chain presented above opens up the possibility of studying the properties of elemen- 
tary excitations of all the general SU(m\n) models in a consistent manner. While we have 
in particular focused on lattice models with inverse square exchange, our results carry over 
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to many other graded permutation models as well. 

This work was supported in part by the U. S. Department of Energy under grant 
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